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Abstract 



We study in detail the barrier distributions extracted from large-angle 
quasi-elastic scattering of heavy ions at energies near the Coulomb barrier. 
Using a closed-form expression for scattering from a single barrier, we compare 
the quasi-elastic barrier distribution with the corresponding test function for 
fusion. We examine the isocentrifugal approximation in coupled-channels cal- 
culations of quasi-elastic scattering and find that for backward angles, it works 
well, justifying the concept of a barrier distribution for scattering processes. 
This method offers an interesting tool for investigating unstable nuclei. We il- 
lustrate this for the ^^Mg + ^''^Pb reaction, where the quadrupole collectivity 
of the neutron-rich ^■^Mg remains to be clarified experimentally. 
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I. INTRODUCTION 



Heavy-ion collisions at energies around the Coulomb barrier provide an ideal opportunity 
to study quantum tunneling phenomena in systems with many degrees of freedom [1,2]. In 
a simple model, a potential barrier for the relative motion between the colliding nuclei is 
created by the strong interplay of the repulsive Coulomb force with the attractive nuclear 
interaction. In the eigenchannel approximation, this barrier is split into a number of dis- 
tributed barriers due to couplings of the relative motion to intrinsic degrees of freedom (such 
as collective inelastic excitations of the colliding nuclei and/or transfer processes), resulting 
in the subbarrier enhancement of fusion cross sections [3]. It is now well known that a bar- 
rier distribution can be extracted experimentally from the fusion excitation function aiy^{E) 
by taking the second derivative of the product Ea^^siE) with respect to the ccnter-of-mass 
energy E, that is, (P{Eai^)/dE'^. This method was first proposed by Rowley, Satchler, 
and Stelson in Ref . [4] , and has stimulated precise measurements of fusion cross sections for 
many systems [5,6] (see Ref. [1] for a detailed review). The extracted fusion barrier distri- 
butions have been found to be very sensitive to the structure of the colliding nuclei, and 
thus the barrier distribution method has opened up the possibility of using the heavy-ion 
fusion reaction as a "quantum tunneling microscope" in order to investigate both the static 
and dynamical properties of atomic nuclei. 

Channel couphngs also affect the scattering process. In Ref. [7], it was suggested that 
the same information as the fusion cross section may be obtained from the cross section 
for quasi-elastic scattering (a sum of elastic, inelastic, and transfer cross sections) at large 
angles. At these backward angles, it is known that the single-barrier elastic cross section 
falls off smoothly from a value close to that for Rutherford scattering at low energies to 
very small values at energies high above the barrier. Timmers et al. therefore proposed to 
use the first derivative of the ratio of the quasi-elastic cross section Uqei to the Rutherford 
cross section ctr with respect to energy, —d{da^c\/ da r) / dE ^ as an alternative representation 
of the barrier distribution [8]. The experimental data of Timmers et al. have revealed [8] 
that the quasi-elastic barrier distribution is indeed similar to that for fusion, although the 
former may be somewhat smeared and thus less sensitive to nuclear structure effects. 

There are certain attractive experimental advantages to measuring the quasi-elastic cross 
section dqci rather than the fusion cross sections dfus to extract a representation of the 
barrier distribution [9]. These are: i) less accuracy is required in the data for taking the 
first derivative rather than the second derivative, ii) whereas measuring the fusion cross 
section requires speciahzed recoil separators (electrostatic defiector/velocity filter) usually 
of low acceptance and efficiency, the measurement of CTqei needs only very simple charged- 
particle detectors, not necessarily possessing good resolution either in energy or in charge, 
and iii) several effective energies can be measured at a single-beam energy, since, in the 
semi-classical approximation, each scattering angle corresponds to scattering at a certain 
angular momentum, and the cross section can be scaled in energy by taking into account 
the centrifugal correction. The last point not only improves the efficiency of the experiment, 
but also allows the use of a cyclotron accelerator where the relatively small energy steps 
required for barrier distribution experiments cannot be obtained from the machine itself. 
This fact was recently exploited by Piasecki et al. [10], who took an astute choice of the 
scattering angles at which (jqei was measured in order to have the energy range necessary. 
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while retaining relatively small energy steps. Moreover, these advantages all point to greater 
ease of measurement with low-intensity exotic beams. 

In this paper, we undertake a detailed discussion of the properties of the quasi-elastic 
barrier distribution. In contrast to the fusion barrier distribution, a theoretical description of 
the quasi-elastic barrier distribution has been limited so far either to a purely classical level 
or to a completely numerical level. Given that many new barrier distribution measurements 
for exotic nuclei are expected to come out in the near future, due to an increasing availability 
of radioactive beams, we believe that it is of considerable importance to clarify the properties 
of the quasi-elastic barrier distribution in a more reliable and transparent way. 

The paper is organized as follows. In Sec. II, we consider a single-barrier system and 
discuss test functions for the barrier distribution, that is the representations of the barrier 
distribution for a single barrier case. We first briefly review the fusion test function, and 
discuss the relation to the barrier penetrability. We then use semi-classical perturbation 
theory [11,12] to derive an analytical expression for the elastic cross section at backward 
angles. Using the formula thus obtained, we discuss the energy dependence of the quasi- 
elastic test function, and compare it with that for the fusion test function. We also discuss 
the scaling property of the quasi-elastic test function obtained at different scattering angles. 
In Sec. Ill, we discuss the barrier distribution for coupled-channels systems. Theoretically, 
barrier distributions have a clear physical meaning only in the limit of zero angular momen- 
tum transfer (that is, in the isocentrifugal approximation [13-21]) with vanishing excitation 
energies for the intrinsic degrees of freedom. In this limit the barrier distribution repre- 
sentation may be derived analytically as a weighted sum of test functions. Nevertheless, a 
simple two-level model suggests that the concept holds to a good approximation even when 
the excitation energy is finite [22]. And of course many experimental data also show well- 
defined barrier structures, which can be reproduced by coupled-channels calculations, even 
for systems where the excitation energies are large. However, although the validity of the 
isocentrifugal approximation has been shown to work well for fusion [18], its applicability 
for scattering processes in the presence of the long-range Coulomb interaction is less clear 
[17,19-21]. We therefore re-examine its validity for the quasi-elastic barrier distribution. In 
Sec. IV, we consider the quasi-elastic barrier distribution as applied to reactions induced by 
exotic nuclei. In particular, we demonstrate its usefulness by showing the possible effects of 
the quadrupole excitation of ^^Mg in the ^^Mg -|- ^osp^ system. We summarize the paper in 
Sec. V. 

II. SINGLE-BARRIER PROBLEMS 

In this section, we discuss heavy-ion reactions between inert nuclei. For such a system, 
the incident fiux of the projectile is either absorbed or elastically scattered from the target 
nucleus. We use a local optical potential which is energy and angular momentum indepen- 
dent. Assuming that the imaginary part of the optical potential is strong enough and is 
localized well inside the Coulomb barrier, the absorption cross section is identified with the 
fusion cross section. 



3 



A. Fusion test function 



Let us first discuss the properties of tlie fusion test function. The classical fusion cross 
section is given by, 

<s{E) ^ nRl {l - ^) e{E-B), (1) 

where Rh and B are the barrier position and the barrier height, respectively. From this 
expression, it is clear that the first derivative of Eaf^^ is proportional to the classical pene- 
trability for a 1-dimensional barrier of height B or eqivalently the s-wave penetrability, 

^[Eatm = nRl e{E - B) ^ nRl P^{E), (2) 
and the second derivative to a delta function. 



dE^ 



[Ea^UE)]=^RiS{E-B). (3) 



In quantum mechanics, the tunneling effect smears the delta function in Eq. (3). An 
analytic formula for the fusion cross section can be obtained if one approximates the Coulomb 
barrier as an inverse parabola, and is given by [23], 



<Tfus(^) = f^R:t In [1 + e'-^^-^^/^^^\ , (4) 

where HQ is the curvature of the Coulomb barrier. Again, the first derivative of Eaf^s is 
proportional to the s-wave penetrability for a parabolic barrier, 

^[Ea,^{E)] = nRl r ! ^ . = nRl P{E). (5) 



dE' ^"^^ ''i + exp[-|g(E-5) 

Defining the function Gfus(-E') as. 



eq. (5) leads to 



GUE) ^ :^^^^[EaUE)l (6) 



r (FY dP{E)_2^ 



where x = 2t:{E — B)/hQ,. This function has the following properties: i) it is symmetric 
around £^ = S, ii) it is centered on E — B, iii) its integral over E is unity, and iv) it has a 
relatively narrow width of around ^r^ln(3 -|- y/S)/n ~ 0.56hn. In the next section, we will 
show that a barrier distribution can be expressed as a weighted sum of normalized functions 
G{E) (see Eq. (21)). The function Gfus(-E), therefore, plays the role of a test function, and 
we call it the fusion test function. 
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B. Quasi-elastic test function 



We now ask ourselves the question of how best to define a similar test function for a 
scattering problem. In the pure classical approach, in the limit of a strong Coulomb field, 
the differential cross sections for elastic scattering at ^ = tt is given by, 

ai{E,ii) = aR{E,i:)e{B-E), (8) 

where aji{E,7r) is the Rutherford cross section. Thus, the ratio af^{E,7r)/aji{E,7r) is the 
classical reflection probability R{E) {— 1 — P{E)), and Eq. (7) suggests that the appropriate 
test function for scattering is [8], 

G,ei(^) = = J • (9) 

In realistic systems, however, due to the effect of nuclear distortion, the differential cross 
section deviates from the Rutherford cross section even at energies below the barrier. Using 
the semi-classical perturbation theory [11,12.24], we derive in the Appendix a semi-classical 
formula for the backward scattering which takes into account the nuclear effect to the leading 
order. The result for a scattering angle 9 reads, 

^^^ = a{E,X^).\S{E,X^)\\ (10) 

where S{E, Ac) is the total (Coulomb -|- nuclear) ^'-matrix at energy E and angular momen- 
tum Ac = rj cot {9/ 2), with rj being the usual Sommerfeld parameter. Note that |5'(£', Ac)|^ 
is nothing but the reflection probability of the Coulomb barrier. For 9 — n, Xc is zero, and 
|5(£;,Ac=0)|2is given by 



exp 




ME-B) 




1 -|- exp 


l-miE-B)] 



\S{E, Ac = 0)1^ = R{E) = ^ (11) 



in the parabolic approximation. a{E, Ac) in Eq. (10) is given by 

V/v(?"c) \/2aTikrj 



a{E, Ac) = 1 + 



ka E 



2VN{rc) - 1 



(12) 



where k = y2fiE/h^, with /i being the reduced mass for the colliding system. The nuclear 

potential Vjv(^c) is evaluated at the Coulomb turning point rc—{r] + ^Jrf -\- X^/k, and a is 
the diffuscncss parameter in the nuclear potential. 

The upper panel of Fig. 1 shows the excitation function of the cross sections at 6* = tt 
for the ^^O -|- ^^^Sm reaction. We use an optical potential of the Woods-Saxon form, with 
parameters Vq — 105.1 MeV, Tq = 1.1 fm, a — 0.75 fm, W=?>Q MeV, rw — 1-0 fm, and aw 
= 0.4 fm. The solid line is the exact solution of the Schrodingcr equation, while the dashed 
line is obtained with the semi-classical formula (10). The dotted line shows the reflection 
probability R{E) = \S{E)\'^. We clearly see that the semi-classical formula accounts well for 
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the deviation of the elastic cross section (7e\{E) from the Rutherford cross section around 
the Coulomb barrier. 

The corresponding quasi-elastic test functions, G^^\{E) = —d/dE{ae\/ (Jr)i are shown in 
the lower panel of Fig. 1. We use a point-difference formula with Ai^cm = 1-8 MeV (as in 
an experiment) in order to evaluate the energy derivative. Notice that the first derivative 
of the reflection probability (dotted line) corresponds to the fusion test function G{^is{E) 
given in Eq. (7). Because of the nuclear distortion factor a(i?, Ac), the quasi-elastic test 
function behaves a little less simply than that for fusion. We find: i) the peak position 
slightly deviates from the barrier height B (by 0.265 MeV for the example shown in Fig. 1), 
and ii) it has a low-energy tail. Eq. (10) indicates that there are two contributions to the 
quasi-elastic test function. One is a{E) ■ dR{E)/dE, and the other da{E)/dE ■ R{E). In 
Fig. 2, we show these two contributions separately. We notice that the low-energy tail of 
the quasi-elastic test function comes from the latter, that is, the energy dependence of the 
nuclear distortion factor a{E). 

Despite these small drawbacks, the quasi-elastic test function Gqei{E) behaves rather 
similarly to the fusion test function Gfus (-£"). In particular, both functions have a similar, 
relatively narrow, width, and their integral over E is unity. We may thus consider that the 
quasi-elastic test function is an excellent analogue of the one for fusion, and we exploit this 
fact in studying barrier structures in heavy-ion scattering. 

Notwithstanding the above comments, it is clear that the quasi-elastic test function 
defined above depends on the scattering angle, and below we shall show how the test function 
can be scaled in terms of an effective energy. 



C. Scaling property of the quasi-elastic test function 

One of the advantages of the quasi-elastic test function over the fusion test function is 
that different scattering angles correspond to the different grazing angular momenta. To 
some extent, the effect of angular momentum can be corrected by shifting the energy by 
an amount equal to the centrifugal potential. Estimating the centrifugal potential at the 
Coulomb turning point Tc, the effective energy may be expressed as [8] 

2//r2 l + sin(e/2) ^ ^ 

In deriving this equation, we have used the definition of Tc, that is, E = ZpZ^e^ + 
A^^^/2/ir^. Therefore, one expects that the function —d/ dE{a^\/ gr) evaluated at an angle 
6 will correspond to the quasi-elastic test function (9) at the effective energy given by eq. 
(13). 

In order to check the scaling property of the quasi-elastic test function with respect to the 
angular momentum. Fig. 3 compares the functions (Jei/cFR (upper panel) and —d/ dE{ae\l (Tr) 
(lower panel) obtained at two different scattering angles. The solid line is evaluated at 6* = tt, 
while the dotted line at = 160°. The dashed line is the same as the dotted line, but shifted 
in energy by E'eff — E. Evidently, the scaling does work well, both at energies below and 
above the Coulomb barrier. 

We should note, however, that as the scattering angle decreases, the scaling becomes less 
good. See Fig. 4 for the scaling property for 9 — 140°. Thus in planning an experiment 
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(especially if it combines data taken in detectors at different angles), one should take careful 
account of this effect. Also at smaller angles, it is well known that the underlying elastic 
cross section will display Fresnel oscillations, which would cause the test function itself (and 
any derived distribution) to oscillate. Detector angles are best chosen to minimise effects of 
Presnel oscillations. 



III. BARRIER DISTRIBUTION FOR MULTI-CHANNEL SYSTEMS 

A. Barrier distributions in the sudden tunneling limit 

Let us now discuss the barrier distributions in the presence of a coupling between the 
relative motion r and an intrinsic degree of freedom ^. The standard way to address the 
effect of the coupling is to solve the coupled-channels equations. For a problem of heavy- 
ion fusion reactions, these equations are often solved in the iso-centrifugal approximation 
[25], where one replaces the angular momentum of the relative motion in each channel 
by the total angular momentum J (this approximation is also referred to as the rotating 
frame approximation or the no-Coriolis approximation in the literature) . The iso-centrifugal 
approximation dramatically simplifies the angular momentum couplings, and reduces the 
dimension of the coupled-channels equations in a considerable way [13-21]. The coupled- 
channels equations in this approximation are given by 

' + \ J + Vo{t) -E + eA ui{t) 



+ E\/^^/W(^™I^Aok7'o)«/'(r) = 0, (14) 



where \(Pim) is an intrinsic wave function which satisfies i^intlv^/M) = ^i\fiM)- We have 
assumed that the coupling Hamiltonian is given by Koup = /('')^A/i(^)^AM(0- The coupled- 
channels equations are solved with the scattering boundary condition for ur{r), 

Mr) - - <; H^j-\kr)Srj, - ^SiHf\k,r) ) , (15) 

where Sj is the nuclear S-matrix. hI \kr) and Hi^\kr) are the incoming and the out- 
going Coulomb wave functions, respectively. The channel wave number ki is given by 
'\j2iji(E — ei)/h^, and ki = k^ — \j2[iElf? . The scattering angular distribution for the 
channel / is then given by [17], 

with 
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where aj{E) and /c(^) are the Coulomb phase shift and the Coulomb scattering amplitude, 
respectively. 

In the limit of e/ 0, the reduced coupled-channels equations (14) are completely 
decoupled. In this limit, the coupling matrix defined as 



Vii' = ei6ij' + y ^ f{r){(pio\Txo\(pi'o) (18) 

can be diagonalized independently of r. It is then easy to prove that the fusion and the 
quasi-elastic cross sections are given as a weighted sum of the cross sections for uncoupled 
eigenchannels [14,15], 

aUE)^j:w^ai:^{E), (19) 

a 

a^,,{E,e) = E^K^) ^Y.w^ai^\E,e), (20) 

I a 

where a^f^^{E) and al,i\E,6) arc the fusion and the elastic cross sections for a potential in 
the eigenchannel a, that is, Va{r) = Vo(r) + Aa(r). Here, Ao(r) is the eigenvalue of the 
coupling matrix (18) (when e/ is zero, AQ,(r) is simply given by • /(r)). The weight factor 
Wa is given by Wa — Uq^, where U is the unitary matrix which diagonalizes Eq. (18). Eqs. 
(19) and (20) immediately lead to the expressions for the barrier distribution in terms of 
the test functions introduced in the previous section, 

DUE) = ^[EaUE)] = E^-'^'R*,. G<'J(E). (21) 

(22) 

As an example of these formulas, let us consider the effect of rotational excitations of 
the target nucleus in the reaction of ^^O with the deformed ^^^Sm. For this problem, cross 
sections (19) and (20) can be computed as [26] 

a(E) ^ d(coseT)(T(E;eT), (23) 

^0 

where 9t is the orientation of the deformed target. The angle dependent potential V{r, 6t) 
is given by, 

V{r, Ot) - Vn{t, Ot) + Vc{t, Ot), (24) 
^"^^ " 1 + exp[(r -R- RtP^Y^o^Ot) - i?r/?4no(^T))] ' ^^^^ 

Voir. ^.) = ^ + E + 'S^l^^^ ^tfr^ (2^) 

Figs. 5(a) and 5(b) show the barrier distributions obtained with Eq. (23) for the fusion 
and the quasi-elastic processes, respectively. We use the potential whose parameters are 
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T/o=220 MeV, i? = 1.1 x {A]^^ + A]^^) fm, and a=0.65 fm. The deformation parameters are 
taken to be /?2=0.306 and /94=0.05. We replace the integral in Eq. (23) with the (/max + 2)- 
point Gauss quadrature [15] with /max=10. That is, we take 6 different orientation angles. 
The contributions from each eigenbarrier are shown by the dashed line in Fig. 5(a) and 
Fig. 5(b). The sohd hne is the sum of all the contributions, which is compared with the 
experimental data [5,8]. The agreement between the calculation and the experimental data 
is reasonable both for the fusion and the quasi-elastic barrier distributions. For the fusion 
barrier distribution -Dfus, the agreement will be further improved if one uses a larger value of 
diffuseness parameter a [5,27] (see the dotted line). Fig. 5(c) compares the fusion with the 
quasi-elastic barrier distributions. These are normalized so that the energy integral between 
50 and 70 MeV is unity. As we discussed in Sec. II for a single barrier case, we see that the 
two barrier distributions show a very similar behavior to each other. 

B. Barrier distributions in systems with finite excitation energy 

In general, the approximation of neglecting the excitation energies ej (that is, the sudden 
tunneling approximation) is valid only for rotational states in heavy deformed nuclei. Despite 
this, however, some of the most interesting effects have been found in the fusion barrier 
distributions for systems involved with highly vibrational nuclei as well [5,6]. One finds that 
the barrier structures still exist, but that the weights of the different barriers can be strongly 
influenced by non-adiabatic effects. In Ref. [22], we have explicitly demonstrated that the 
fusion cross sections are in general given by Eq. (19), but with the energy dependent 
weight factors Wa{E) (in the sudden tunneling limit, the weight factors become energy 
independent). For a simple two-channel problem, we found that although the weights may 
depend strongly on the excitation energy, their dependence on the incident energy is weak, 
suggesting that the concept of a barrier distribution holds good even for flnite intrinsic 
excitation energies [22]. Since the quasi-elastic barrier distribution Gqei{E) is related to the 
fusion barrier distribution Gius{E) through flux conservation (unitarity), a similar situation 
can be expected for the quasi-elastic barrier distribution. 

C. Applicability of the iso-centrifugal approximation 

As we have mentioned in Sec. I, the validity of the iso-centrifugal approximation has 
been well tested for heavy-ion fusion reactions [18]. In contrast, it is known that the approxi- 
mation fails to reproduce the exact result for scattering angular distributions in the presence 
of the long-range Coulomb force. The effect of the coupling is somewhat overestimated in 
the isocentrifugal approximation, and simple recipes to renormalize the coupling strength 
have been proposed in order to cure this problem [17,19-21]. On the other hand, Esbensen 
et al. have argued, based on semi-classical considerations, that the isocentrifugal approxi- 
mation (without renormalization of the coupling strength) works better for backward angle 
scattering [17]. 

Since it has not yet been clear how well the isocentrifugal approximation works in con- 
nection with the quasi-elastic barrier distribution, we re-examine in this subsection the 
performance of the approximation for large-angle scattering. To this end, we consider the 
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effect of quadrupole plionon excitations in tfie target nucleus for the ^^O + ^ Sm reaction. 
In order to emphasise the couphng effect, we increase the couphng strength and reduce the 
excitation energy from the physical values. The values which we use are: (32=0.2 (with rcoup 
= 1.06 fm) and 62=0.5 MeV. We have checked that our conclusions are not altered irrespec- 
tive of the values of P2 and €2- For simplicity, we consider only a single phonon excitation, 
and employ the linear coupling approximation [28] . We use the same optical potential as in 
Sec. 11. 

Figure 6 shows the partial cross sections at £^cm=65 MeV for the angle-integrated inelastic 
scattering (upper panel) and for the fusion reaction (lower panel) as a function of the initial 
orbital angular momentum = J. The solid line is the exact result of the coupled-channels 
equations with the full angular momentum couphngs, while the dashed line is obtained with 
the iso-centrifugal approximation. We find that the isocentrifugal approximation works 
rather well for J < 20, although the agreement is poor for larger values of J. For fusion, 
only small values of J contribute, and the isocentrifugal approximation always makes an 
excellent approximation. Fig. 7 shows the angular distributions for the elastic (upper 
panel) and inelastic scattering (lower panel). Although the isocentrifugal approximation 
does not reproduce the main structure of the angular distribution, it indeed works very 
nicely at backward angles where the main contribution comes from small values of angular 
momentum (see Eq. (10) and Fig. 6). In fact, the isocentrifugal approximation almost 
reproduces the exact result for the scattering angles 9cm > 130°. 

Fig. 8 shows the excitation function for quasi-elastic scattering (upper panel) and its 
energy derivative calculated at = 170° in the laboratory frame. One sees that the isocen- 
trifugal approximation well reproduces the exact solution. We thus conclude that the the 
isocentrifugal approximation works sufficiently well for studies of quasi-elastic barrier dis- 
tributions. This fact not only makes the coupled-channels calculations considerably easier, 
but also assures the similarity of fusion and quasi-elastic distributions even in the presence 
of channel couplings. 

IV. QUASI-ELASTIC SCATTERING WITH RADIOACTIVE BEAMS 

It has been well recognized that low-energy reactions provide an ideal tool to probe 
the detailed structure of atomic nuclei. The heavy-ion fusion reaction around the Coulomb 
barrier is one of the typical examples. In the last decade, many high-precision measurements 
of fusion cross sections have been made, and the nuclear structure information has been 
successfully extracted through the representation of the fusion barrier distribution [1]. 

Low-energy radioactive beams have also become increasingly available in recent years, 
and heavy-ion fusion reactions involving neutron-rich nuclei have been performed for a few 
systems [29-33]. New generation facilities have been under construction at several labora- 
tories, and many more reaction measurements with exotic beams at low energies will be 
performed in the near future (see Ref. [34] for a recent theoretical review). Although it 
would still be difficult to perform high-precision measurements of fusion cross sections with 
radioactive beams, the measurement of the quasi-elastic barrier distribution, which can be 
obtained much more easily than the fusion counterpart as we mentioned in the introduction, 
may be feasible. Since the quasi-elastic barrier distribution contains similar information as 
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the fusion barrier distribution, the quasi-elastic measurements at backward angles may open 
up a novel way to probe the structure of exotic neutron-rich nuclei. 

In order to demonstrate the usefulness of the study of the quasi-elastic barrier distribution 
with radioactive beams, we take as an example the reaction ^^Mg and ^°^Pb. The neutron- 
rich ^^Mg nucleus has attracted much interest as evidence for the breaking of the A^=20 
spherical shell closure. In this nucleus, a large B[E2) value (454±78 e^fm'^ [35] and 622± 
90 e^fm'^ [36]) and a small value of the excitation energy of the first 2+ state (885 keV) 
[35] have been experimentally observed. The authors of Refs. [35-37] argue that these large 
collectivities may be indicative of a static deformation of ^^Mg. On the other hand, mean- 
field calculations [38] as well as quasiparticle random-phase approximation (QRPA) [39] with 
the Skyrme interaction suggest that "^^Mg may be spherical. In fact the energy ratio between 
the first 4+ and the first 2+ states, E^+jE^+i is 2.6 [37], which is between the vibrational 
and rotational limits [39]. 

Note that the distorted-wave Born approximation (DWBA) yields identical results for 
both rotational and vibrational couplings (to first order). In order to discriminate whether 
the transitions are vibration-like or rotation-like, at least second-step processes (reorienta- 
tion and/or couplings to higher members) are necessary. The coupling effect plays a more 
important role in low-energy reactions than at high and intermediate energies. Therefore, 
quasi-elastic scattering around the Coulomb barrier may provide a useful method of clarify- 
ing the nature of the quadrupole collectivity of ^^Mg. 

Fig. 9 shows the excitation function of the quasi-elastic scattering (upper panel) and the 
quasi-elastic barrier distribution (lower panel) for this system. The solid and dashed lines 
are results of coupled- channels calculations where ^^Mg is assumed to be a rotational or a 
vibrational nucleus, respectively. We estimate the coupling strength from the measured 
B{E2) value [35] to be 0.51. We include the quadrupole excitations in ^^Mg up to the 
second member (that is, the first 4+ state in the rotational band for the rotational coupling, 
or the double phonon state for the vibrational coupling). In addition, we include the single 
octupole phonon excitation at 2.615 MeV in ^°^Pb [40]. The potential parameters which 
we use are Vo=180 MeV, ro=1.15 fm, and a=0.63 fm, that give the same barrier height 
(5=106.6 MeV) as the Akyiiz-Winther potential [41]. For the imaginary potential, we use 
1^=50 MeV, Tyj—l.^ fm, and a^=0.4 fm, but the results are insensitive to this as long as 
it is localized inside the barrier with a large enough strength. We use the computer code 
CQUEL [42] in order to integrate the coupled-channels equations. This code is a version of 
CCFULL [25] , where the coupling is treated to all orders in the coupling hamiltonian and the 
isocentrifugal approximation is employed in order to reduce the dimension of the coupled- 
channels equations. In the code CQUEL, we use the regular boundary condition at the origin, 
instead of the incoming boundary condition, and we remove the restriction of CCFULL, which 
computes only the fusion cross sections. 

In the figure, we can see well separated peaks in the quasi-elastic barrier distribution both 
for the rotational and for the vibrational couplings. Moreover, the two fines are consider- 
ably different at energies around and above the Coulomb barrier, although the two results 
are rather similar below the barrier. We can thus expect that the quasi-elastic barrier dis- 
tribution can indeed be utilized to discriminate between the rotational and the vibrational 
nature of the quadrupole collectivity in ^^Mg, although these results might be somewhat per- 
turbed by other effects which are not considered in the present calculations, such as double 
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octupole-phonon excitations in the target, transfer processes or hexadecapole deformations. 



V. SUMMARY 



The quasi-elastic barrier distribution is a counterpart of the fusion barrier distribution 
in the sense that the former is related to the reflection probability of a potential barrier 
while the latter is related to the transmission. In this paper, we have studied some detailed 
properties of the quasi-elastic barrier distribution. Using semi-classical perturbation theory, 
we have obtained an analytic formula for the quasi-elastic barrier distribution for a single 
barrier (that is, the quasi-clastic test function). The formula indicates that this test function 
consists of two factors: one is related to the effect of the nuclear distortion of the classical 
trajectory, while the other is the reflection probability of the potential barrier. Due to the 
nuclear distortion, we found that the quasi-elastic barrier distribution is slightly less well 
behaved than the fusion barrier distribution. For instance, the peak position of the quasi- 
elastic barrier distribution slightly deviates from the barrier height, and it has a low-energy 
tail. Nevertheless, the quasi-elastic barrier distribution behaves rather similarly to that for 
fusion on the whole, and both are sensitive to the same nuclear structure effects. 

In mult i- channels systems, the validity of the barrier distribution relies on the isocentrifu- 
gal approximation, where the angular momentum of the relative motion in each channel is 
replaced by the total angular momentum J. We have examined the applicability of this ap- 
proximation for scattering processes and have found that it works well at least for backward 
angles, where such experiments are performed. 

The measurement of quasi-elastic barrier distributions is well suited to future experiments 
with low-intensity exotic beams. To illustrate this fact, we have discussed as an example, the 
effect of quadrupole excitations in the neutron-rich ^^Mg nucleus on quasi-elastic scattering 
around the Coulomb barrier, and argued that the quasi-elastic barrier distribution would 
provide a useful tool to clarify whether ^^Mg is spherical or deformed. In this way, we expect 
that the barrier distribution method will open up a novel means to allow the detailed study 
of the structure of neutron-rich nuclei in the near future. 



We thank E. Piasecki and E. Crema for helpful discussions. 

APPENDIX A: SEMI-CLASSCIAL PERTURBATION THEORY 

In this Appendix, we derive Eq. (10) for the backward- angle elastic cross section using 
semi-classical perturbation theory. Our formula is an improvement of the one in Ref. [12], 
since we take into account the effect of nuclear distortion of the classical trajectory [24]. 

The scattering amplitude f{9) for a spherical optical potential is given by. 
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where 9 is the scattering angle and k — Since we are interested in backward 

scattering near ^ ~ tt, we replace the Legendre polynomials P;(cos^) with their asymptotic 
form, 



sin^ 



(A2) 



where Jo{0) is the Bessel function. We now apply the well known Poisson sum formula to 
Eq. (Al) to obtain 



m = T 



I h-e 



k\ sin 9 



E( 



d\ \S{\)Jo{\{n-9))e 



(2n-l)i7rA 



(A3) 



where X — I + 1/2. At energies around the Coulomb barrier and for backward scattering, 
the contribution from n = dominates the sum in Eq. (A3) [11]. Taking only n = and 
evaluating the integral in the stationary phase approximation, one obtains (see Sec. 5.7 of 
Ref. [11]), 



A 



A;2sin^|e'(A)| 



(A4) 



where 0(A) = 2Re6'{X) is the deflection function, 6{X) being the phase shift, and A satisfies 
the stationary phase condition ©(A) = 9. Here, the dash denotes the derivative with respect 
to the argument. This equation yields 



a{9l 
(Tr{9) 



e'(A) 



(A5) 



Landowne and Wolter evaluated Eq. (A5) using a perturbation theory based on the 
semiclassical approximation [12]. The stationary condition 0(A) = Oc(Ac) = 9 and the 
definition of the nuclear defiection function, 0(A) = 0c(A) + 0Ar(A), yield [12] 



0'c(Ac 



0'(A) 



£^i + _^e^(A,) + |e'^(Ae), 



(A6) 



to first order in A — Ac. In deriving this equation, we have assumed that rj is much larger 
than Ac. In the semiclassical approximation, the nuclear phase shift is given by [11] 



/•oo /"OO 

Sn{X) = / k{r)dr — / kc{r)dr, 

Jri Jtc 

kir) = ^J2^iiE - T/jv(r) - Vdr) - Vx{r))/h\ 

Ur) = ^^2^^{E~Vc{r)-Vx{r))/n\ 



(AT) 

(A8) 

(A9) 



where Vn{t) and Vc{t) are the nuclear and the Coulomb potentials, respectively, and V\{r) = 
A^/i^/2/ir^ is the centrifugal potential. The classical turning points ri and Tc satisfy k(ri) = 
kci^c) = 0. To first order in the nuclear potential, the semi-classical phase shift is given by 
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Expanding kc{r) around r — rc and assuming that VN{r) ~ — Voe~^/" near Tc, one obtains 
[11,12] 

25^{\) -y^(r,)^^^ + 0{Xyv'). (All) 

Using the perturbativc phase shift (All) in Eq. (A6), Landowne and Wolter obtained a 
simple form for the backward cross sections which is given by [12], 

An improved formula may be obtained by taking into account the effect of nuclear dis- 
tortion of the classical trajectory. To this end, we follow the method suggested by Brink 
and Satchler [24]. Transforming the coordinate in the first integral in Eq. (A7) to the one 
which satisfies k{r) = kc{s), the semi-classical phase shift may be expressed as 

fOO rj fOO rj 

5n{\)= kc{s)-{r{s) - s)ds ^ - {r{s)-s)-kc{s)ds. (A13) 
Jrc (Is Jtc (Is 

The condition k{r) = kc{s) yields = Vn{s) + {V^{s) + ^^(s) + Vi{s)){r - s) to first order 
in r — s. We thus obtain 

^"^^^ ^ -¥ L vM + vs{s) + viis) ■ M?''' ^ ^ 



V6{rc) + Viivc) 



n roo -VN(r) , / A 

■ / , J ; dr. (A15) 

h^Jr. kjr) ^ ' 



-Vjv(^c) \ V2a7ikr] 



E 



+ 0{\^/if). (A16) 



Here, we have expanded r — s with respect to V^r in Eq. (A14) and evaluated it at the radius 
rc- Substituting Eq. (A16) into Eq. (A5), we finally obtain Eq. (10). 

Fig. 10 compares the semi-classical formula with the exact result (solid line) for the 
^^O + ^"^"^Sm reaction. We use the same optical potential as in Sec. 11. The dotted line is 
obtained by the semi-classical perturbation of Landowne and Wolter, Eqs. (All) and (A12). 
The dashed line is the result of semi-classical approximation which takes into account the 
nuclear distortion, Eqs. (A16) and (10). We see that the semi-classical perturbation theory 
works reasonably well around the Coulomb barrier when the effect of nuclear distortion is 
included. The deviation of the nuclear phase shift from the exact solution above the barrier 
would be improved by using the full semi-classical phase shift [43]. However, we note that 
the backward cross sections are already reproduced reasonably well even with the present 
semi-classical perturbation theory. 
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FIG. 1. The ratio of elastic scattering to the Rutherford cross section &t 6 = tt (upper panel) 
and the quasi-elastic test function Gqei(-E) = —d/dE{aei/crji) (lower panel) for the ^^O + ^^^Sm 
reaction. The solid line is the exact solution of the optical potential, while the dashed line is 
obtained with the semi-classical perturbation theory. The dotted line denotes the reflection prob- 
ability R{E) = \S{E)\^ for s-wave scattering. 
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FIG. 2. Two separate contributions to the quasi-elastic test function. The sohd hne shows the 
function a{E) ■ dR{E)/dE, while the dashed line shows da{E)/dE ■ R{E), where a{E) and R{E) 
are the nuclear distortion function and the reflection probability, respectively. 
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FIG. 3. Comparison of the ratio dcxja^ (upper panel) and its energy derivative —d/ dE{ac\/ (Tr) 
(lower panel) evaluated at two different scattering angles. The solid line is for 9 = n, while the 
dotted line is for 6 = 160". The dashed line is the same as the dotted line, but is shifted in energy 
by an amount equal to the centrifugal potential evaluated at the distance of closest approach of 
the Rutherford trajectory. 
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FIG. 4. The same as fig.3, but for = 140°. 
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FIG. 5. (a) The fusion barrier distribution £>fus(-E) = (f{Eat^)/dE'^ for the ^^^O + ^^^^Sm 
reaction. The soUd Une is obtained with the orientation-integrated formula with (32 = 0.306 and 
P4= 0.05. The dashed lines indicate the contributions from the six individual eigenbarriers. These 
lines are obtained by using a Woods-Saxon potential with a surface diffuseness parameter a of 0.65 
fm. The dotted line is the fusion barrier distribution calculated with a potential which has a = 
1.05 fm. Experimental data are taken from Ref. [5]. (b) Same as Fig. 5(a), but for the quasi-elastic 
barrier distribution Dqei{E) = —d[aqe\{E,TT)/aji{E,Tr)]/dE. Experimental data are from Ref. [8]. 
(c) Comparison between the barrier distribution for fusion (solid line) and that for quasi-elastic 
scattering (dashed line). These functions are both normalized to unit area in the energy interval 
between 50 and 70 MeV. 
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FIG. 6. The effect of a quadrupole-phonon excitation in the target nucleus on the partial cross 
sections for the ^^O + ^^^Sm reaction at £'cm=65 MeV. The upper and the lower panels show the 
angle-integrated inelastic scattering and the fusion cross sections, respectively. The solid line is 
the solution of the coupled-channels equations with the full angular momentum coupling, while the 
dashed line is obtained in the isocentrifugal approximation. 
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FIG. 7. The angular distributions for the elastic (upper panel) and the inelastic (lower panel) 
scattering for the ^^O + ^^^Sm reaction at £^cm=65 MeV. The significance of each line is the same 
as in Fig. 5. 
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FIG. 8. The excitation function for quasi-elastic scattering (upper panel) and the quasi-elastic 
barrier distribution (lower panel) for the ^^O -|- ^^^Sm reaction calculated dX 9 = 170° in the 
laboratory frame. The significance of each line is the same as in Fig. 5. 



24 



0.8- 



^ 0.6 



D 0.4 



0.2 h 





32-. 208_, 

Mg + Pb 



Rotational 
Vibrational 




95 100 105 110 115 120 125 



E (MeV) 

cm ^ ^ 

FIG. 9. The excitation function for quasi-elastic scattering (upper panel) and the quasi-elastic 
barrier distribution (lower panel) for the ^^Mg + ^°^Pb reaction around the Coulomb barrier. The 
solid and the dashed lines are the results of coupled-channels calculations which assume that 22 Mg 
is a rotational and a vibrational nucleus, respectively. The single octupole-phonon excitation in 
208p|-, jg g^jgQ included in the calculations. 
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FIG. 10. Comparison of the semi-classical formulae with the exact solution for the ^''O+^'^^Sm 
reaction. The upper and the lower panels show the nuclear phase shift and the ratio of the elastic 
to the Rutherford cross sections at the scattering angle tt, respectively. The solid line is obtained by 
numerically integrating the Schrodinger equation, while the dotted line is the result of the primitive 
semi-classical perturbation theory, Eqs. (All) and (A12). The dashed line indicates the result of 
the semi-classical perturbation theory which takes into account the effect of nuclear distortion of 
the classical trajectory, Eqs. (A16) and (10). 
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